Elastic Properties of C and B^CLN^ composite nanotubes 
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We present a comparative study of the energetic, structural and elastic properties of carbon and 
composite single-wall nanotubes, including BN, BCs and BC2N nanotubes, using a non-orthogonal 
tight-binding formalism. Our calculations predict that carbon nanotubes have a higher Young 
modulus than any of the studied composite nanotubes, and of the same order as that found for 
defect-free graphene sheets. We obtain good agreement with the available experimental results. 



Carbon nanotubes [GJ were first discovered by Iijima 
in the early nineties, as a by-product of fullerene synthe- 
sis. Since then there has been an ever-increasing interest 
in these new forms of carbon, partly due to their novel 
structures and properties, but perhaps more so due to 
the wealth of potentially important applications in which 
nanotubes could be used. Indeed many applications have 
already been reported, from their use as atomic-force mi- 
croscope tips [|| , to field emitters \ 4 ] , nanoscale electronic 
devices || or hydrogen storage [6|, to cite a few. But 
probably the highest potential of nanotubes is in connec- 
tion with their exceptional mechanical properties. 

After the discovery of graphitic nanotubes it was pos- 
tulated that other compounds forming laminar graphite- 
like structures could also form nanotubes Jul. In partic- 
ular, BN g, BC 3 @, BC 2 N 0] and CN (n) nanotubes 
were predicted on the basis of theoretical calculations. 
BN @, BC 3 and BC 2 N Q have since been synthe- 
sized, though some uncertainty remains as to the actual 
structure of BC2N nanotubes |l4| ]. 

The mechanical properties of carbon nanotubes have 
been the subject of a number of theoretical p5|-pl[| as 
well as experimental p2|-|23|] studies. On the theoreti- 
cal side, studies have been mostly carried out using em- 
pirical potentials, although Molina et al. jl7) employed 
an orthogonal tight-binding model in their work. The 
most extensive theoretical study of the elastic proper- 
ties of carbon nanotubes to date is that of Lu |2l[] , who 
used an empirical pair potential model to estimate the 
Young modulus, Poisson ratio and other elastic constants 
of both single-wall and multi-wall nanotubes, as well as 
nanotube ropes. However, it was not possible to ex- 
tend this study to composite nanotubes, given that no 
empirical potential models akin to that used for carbon 
exist for the composite systems. The behavior of car- 
bon nanotubes subject to large axial strains has been 
studied by Yakobson et al. |l8). The bending of car- 
bon nanotubes has been studied experimentally and us- 
ing simulation techniques by Iijima et al. ppj ] . The Young 
modulus of carbon multi-wall nanotubes has been exper- 
imentally determined by Treacy et al. [ p2| using thermal 
vibration analysis of cantilevered tubes. They obtained 
a mean value of the Young modulus of 1.8 ± 1.4 TPa. 



More recently, Wong et al. 23 have obtained a value of 
1.28±0.59 TPa, by recording the force needed to bend an- 
chored nanotubes using an AFM. Chopra and Zettl |24|] 
have also used thermal vibration analysis to estimate the 
Young modulus of multi-wall BN nanotubes, obtaining a 
value of 1.22 ± 0.24 TPa. These experimental and the- 
oretical studies confirm the expectation that nanotubes 
have exceptional stiffness, and could therefore be used in 
the synthesis of highly resistant composite materials. 

In this work we study the structural, energetic and 
mechanical properties of both carbon and composite nan- 
otubes, paying special attention to the mechanical prop- 
erties, since these are expected to play such an impor- 
tant role in many practical applications. This is the 
first time that such a detailed comparative study is un- 
dertaken. In the majority of the calculations reported 
here the atomic interactions have been modeled using a 
non-orthogonal tight-binding scheme due to Porezag and 
coworkers [^5|. Tight-binding (TB) methods |2(| offer a 
good compromise between the more accurate but much 
more costly first principles |27j] techniques, and empirical 
potentials p8j, which are cheaper to use, but often not 
transferable to configurations different to those for which 
they have been fitted. 

In the TB scheme used here, the hopping integrals 
used to construct the Hamiltonian and overlap matri- 
ces are tabulated as a function of the internuclear dis- 
tance on the basis of first-principles density-functional 
theory (DFT) calculations employing localized basis sets, 
retaining only one- and two-center contributions to the 
Hamiltonian matrix elements |2f|. A minimal basis set 
corresponding to a single atomic-like orbital per atomic 
valence state is used. More details on the construction of 
the TB parametrisation used in this work can be found 
in ref. §§. 

Using the non-orthogonal TB scheme briefly outlined 
above we have performed a series of calculations aimed 
at characterizing the elastic properties of single- wall nan- 
otubes. In particular, we have considered C, BN, BC3 
and BC2N (n,n) and (n,0) (i.e. non-chiral) nanotubes. 
Two structures having the same stoichiometry are possi- 
ble for the BC 2 N nanotubes. Only the structure known 
as II p0|,p9|] is considered here, as this is predicted to 
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be the most stable of the two. In addition, we have per- 
formed calculations for the chiral (10,5) and (10,7) C nan- 
otubes. We have also carried out Plane Wave (PW) pseu- 
dopotential DFT calculations within the local density ap- 
proximation (LDA) for the (6,6) C and BN nanotubes, 
for comparison with the TB results. Our PW calcula- 
tions used Troullier-Martins pseudopotentials |3(]] . A cut 
off of 40 Ry was used in the basis set, and 10 Monkhorst- 
Pack [^7j points to sample the one-dimensional Brillouin 
zone. The hexagonal supercell was chosen large enough 
so as to ensure that the minimum distance between a 
tube and any of its periodic images was larger than 5.5 A. 
Our TB calculations were performed using T-point sam- 
pling only, but using periodically repeated cells which 
where large enough along the axial direction so as to en- 
sure that total energy differences were converged to an 
accuracy approximately equal to that achieved with the 
PW calculations. 

The Poisson ratio is defined via the equation 
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Here, e is the axial strain, R eq is the equilibrium tube 
radius, R is the tube radius at strain e and a is Poisson's 
ratio. The values of a obtained for a number of repre- 
sentative tubes considered in this work are reported in 
Table |. Regarding the Young modulus, its conventional 
definition is 
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where Vq is the equilibrium volume, and E is the strain 
energy. In the case of a single- wall nanotube, this defini- 
tion requires adopting a convention in order to define Vo, 
which for a hollow cylinder is given by Vo = 2irLR5R, 
where L is the length, R the radius and 5R is the shell 
thickness. Different conventions have been adopted in the 
past; for example, Lu [^l| recently took SR — 0.34 nm, 
i.e. the interlayer separation in graphite, while Yakobson 
et al Q took the value SR = 0.066 nm. We follow a dif- 
ferent path. Rather than adopting an ad hoc convention, 
we use a different magnitude to characterize the stiffness 
of a single-wall nanotube, which is independent of any 
shell thickness. We define 
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Here, S$ is the surface defined by the tube at equilibrium. 
The value of the Young modulus for a given convention 
value SR is given by Y = Y S /5R. In Table Q we give the 
values obtained for Y, for a number of tubes. 
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FIG. 1. Curvature strain energy as a function of the equi- 
librium tube diameter, as obtained from the tight-binding cal- 
culations, for C, BN and BC3 nanotubes. ab initio results are 
from refs. [8-10]. 

Our first observation is that PW and TB results give 
very similar answers for all the calculated properties. The 
differences in values of the Young modulus (calculated 
adopting the convention SR = 0.34 nm) for the (6,6) C 
nanotube are of the order of 0.12 TPa. 
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(n,m) 


D eg (nm) 


a 


Y s (TPa • nm) Y (TPa) 


C 


(10,0) 


0.791 


0.275 


0.416 


1.22 




(6,6) 


0.820 


0.247 


0.415 


1.22 






(0.817) 




(0.371) 


(1.09) 




(10,5) 


1.034 


0.265 


0.426 


1.25 




(10,7) 


1.165 


0.266 


0.422 


1.24 




(10,10) 


1.360 


0.256 


0.423 


1.24 




(20,0) 


1.571 


0.270 


0.430 


1.26 




(15,15) 


2.034 


0.256 


0.425 


1.25 


BN 


(10,0) 


0.811 


0.232 


0.284 


0.837 




(6,6) 


0.838 


0.268 


0.296 


0.870 






(0.823) 




(0.267) 


(0.784) 




(15,0) 


1.206 


0.246 


0.298 


0.876 




(10,10) 


1.390 


0.263 


0.306 


0.901 




(20,0) 


1.604 


0.254 


0.301 


0.884 




(15,15) 


2.081 


0.263 


0.310 


0.912 


BC 3 


(5,0) 


0.818 


0.301 


0.308 


0.906 




(3,3) 


0.850 


0.289 


0.311 


0.914 




(10,0) 


1.630 


0.282 


0.313 


0.922 




(6,6) 


1.694 


0.279 


0.315 


0.925 


BC 2 N II 


(7,0) 


1.111 


0.289 


0.336 


0.988 




(5,5) 


1.370 


0.287 


0.343 


1.008 



TABLE I. Structural and elastic properties of selected nan- 
otubes obtained from the tight-binding calculations reported 
here. Young modulus values given in parenthesis were ob- 
tained from first-principles calculations. Also the value of Y 
with the convention SR = 0.34 nm is given for comparison. 
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It is worth pointing out that this difference is small 
compared to the uncertainty with which Y can at 
present be experimentally determined for multi-wall nan- 
otubes p^-p^]. Structural properties deduced from the 
PW and TB calculations are also in very good agreement 
for both C and BN systems. The equilibrium diameter, 
as seen m Table |, differs by about 1 % or less. The near- 
est neighbor distance in the case of the (6,6) C nanotube 
is 1.42 A in both the PW and TB calculations. For the 
BN (6,6) nanotube, the results are 1.43 A and 1.45 A for 
the PW and TB calculations, respectively. 

An interesting magnitude associated with nanotubes is 
the curvature energy or strain energy E s , which we define 
as the difference of the energy per atom in the tube and 
that in the corresponding infinite flat sheet. In Fig. we 
plot the strain energy obtained from our TB calculations 
for C, BN and BC3 nanotubes as a function of the tube 
diameter. It can be seen that the characteristic behavior 
E s oc l/D 2 , where D is the tube diameter, is obtained. 
Our calculations indicate that the strain energy at a given 
tube diameter is highest for C nanotubes, and that both 
BN and BC3 nanotubes have very nearly the same strain 
energy. The fact that these composite nanotubes have 
smaller strain energy than pure carbon nanotubes is in 
agreement with previous first-principles results MM. 
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FIG. 2. Buckling in the BN nanotube equilibrium struc- 
tures vs. tube diameter. We define the buckling as the mean 
radius of the nitrogen atoms minus the mean radius of the 
boron atoms. The ab initio results are from ref. [8]. 

A structural feature which is specific to the BN nan- 
otubes is the presence of a certain degree of buckling on 
the tube surface, which results from the B atoms displac- 
ing towards the tube axis, while the N atoms displace in 
the opposite direction. Notice the good agreement with 
ab initio results H . As for the strain energy, the buckling 
effect decreases rapidly with increasing tube diameter, 
going to the flat BN sheet limit of zero buckling. This 
tendency of BN nanotubes to buckle, which is a result 



of the slightly different hybridizations of B and N in the 
curved hexagonal layer, will have the effect of forming a 
surface dipole, a fact that could be relevant for potential 
applications of these tubes. 

In Fig. H we have plotted the values of Y s obtained for 
the different tubes. The first feature to be noticed is the 
fact that both for (n,n) and (n,0) nanotubes, the carbon 
tubes are predicted to be significantly stiffer than any 
of the composite tubes. The BC2N are predicted to be 
somewhat stiffer than BN and BC3 tubes. The value of 
0.43 TPa-nm obtained for the widest C nanotubes corre- 
sponds to a Young modulus of 1.26 TPa, taking a value 
of 0.34 nm for the graphene sheet thickness. This value 
is in excellent agreement with the experimental result of 
1.28 ±0.59 TPa of Wong et al. ||. Although our results 
are only for single- wall tubes, it can be expected that 
the elastic properties of multi-wall tubes and nanoropes 
be mostly determined by the strength of the C-C bonds 
in the bent graphene sheets, and thus be very similar to 
those of single- wall tubes. Our results for C nanotubes 
are also in reasonable agreement with the measurements 
of Treacy et al. @ (1.8±1.4 TPa). Chopra and Zettl ||] 
obtain a value of 1.22 TPa with an estimated 20 % error 
for multi-wall BN nanotubes. This value is again some- 
what larger than what we obtain for BN nanotubes, but 
nevertheless, the agreement is close. Lu's |^l[ estimation 
of the Young modulus for single-wall C nanotubes gives 
results which are slightly smaller than ours (0.97 TPa), a 
difference which is most likely due to the different models 
used in his calculations and ours. 

Our calculations predict that there is a small depen- 
dence of Y s on the tube diameter, but this dependence 
is noticeable only for small values of the tube diameter, 
the limiting (diameter independent) value being rapidly 
obtained at the range of experimentally observed single- 
wall tube diameters (~ 1.2 nm). This is in contrast to 
the results of Lu [^lj , which are almost completely inde- 
pendent of the tube size and chirality. We believe that 
the apparent insensibility of the Young modulus on the 
tube size and chirality observed by Lu is due to the fact 
that an empirical pair potential was used in his calcula- 
tions, and such a model will not reflect the effects that 
the curvature will have on the bonding properties of the 
system. In the limit of large tube diameters, we could 
expect that the elastic properties would correspond to 
those of a plane, defect-free, graphitic sheet. Indeed, cal- 
culations of Y s for plane graphene and BN sheets give 
0.41 and 0.30 TPa • nm respectively, which can be seen 
to be very similar to the results obtained for C and BN 
nanotubes of the largest diameter we studied. It is worth 
noticing that the limiting value of the Young modulus 
as a function of tube diameter is reached from below, 
which is consistent with the expectation that tubes of 
higher curvature (i.e. smaller diameter) will have weaker 
bonds, which would result in a slight reduction of the 
Young modulus. 
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FIG. 3. Young modulus as a function of the tube diameter 
for C, BN, BC3 and BC2N (structure II only), as calculated 
from the tight-binding simulations. Results obtained for (n,n) 
nanotubes (filled symbols) , (n,0) nanotubes (empty symbols) 
and also for C (10,5) (+) and (10,7) (x) are shown. 

To summarize, we have carried out an extensive study 
of the energetic, structural and elastic properties of 
both graphitic and composite nanotubes, using a non- 
orthogonal Tight-Binding scheme. The agreement ob- 
tained between the TB results and the first-principles 
calculations reassure us in our conclusion that the TB 
model employed here gives a good description of the stud- 
ied features of nanotubes. We have obtained good agree- 
ment with the existing experimental measurements of the 
Young modulus for multi-wall C and BN nanotubes. Our 
results indicate that graphitic nanotubes are stiffer than 
any of the composite nanotubes considered in this work, 
and that the elastic properties of single-wall nanotubes 
are of the same order of magnitude as those of the cor- 
responding flat sheets. Although the BN nanotubes are 
predicted to have a somewhat smaller Young modulus 
than the C nanotubes, they remain considerably stiff. 
This fact, combined with their insulator character ||] 
makes them suitable for applications in which electrically 
insulating high-strength materials are needed. 
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